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FOREWORD 


This  report  is  one  of  four  reports  to  be  prepared  by  Structural 
Mechanics  Associates  under  Navy  Contract  No.  N156-46654.  This  con¬ 
tract  was  initiated  under  Work  Unit  No.  530/07,  "Development  of 
Optimization  Methods  for  the  Design  of  Composite  Structures  Made 
from  Anisotropic  Material"  (1-23-96)  and  was  administered  under  the 
direction  of  the  Aeronautical  Structures  Laboratory,  Naval  Air  Engi¬ 
neering  Center,  with  Messrs.  R.  Molella  and  A.  Manno  acting  as  Pro¬ 
ject  Engineers.  The  reports  resulting  from  this  contract  will  be 
forwarded  separately.  Three  reports  are  completed  and  cover  work 
from  '4  May  1965  to  31  December  1966.  The  title  and  approximate 
forwarding  date  for  each  report  are  as  follows: 

NAEC-ASL-1109,  "Structural  Optimization  of  Corrugated  Core 
and  Web  Core  Sandwich  Panels  Subjected  to  Uniaxial  Compression," 
dated  15  May  1967.  Forwarding  date,  June  1967. 

NAEC-ASL-1110,  "Structural  Optimization  of  Flat,  Corrugated 
Core  and  Web  Core  Sandwich  Panels  Under  In-Plane  Shear  Loads 
and  Combined  Uniaxial  Compression  and  In-Plane  Shear  Loads," 
dated  1  June  1967.  Forwarding  date,  July  1967. 

NAEC-ASL-1111,  "A  Method  for  Weight  Optimization  of  Flat  Truss 
Core  Sandwich  Panels  Under  Lateral  Loads,"  dated  15  June  1967. 
Forwarding  date,  July  1967. 
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SUMMARY 

This  report  presents  a  method  for  optimizing  on  &  weight  basis  flat 
truss  core  sandwich  panels  under  lateral  loads. 

To  solve  this  type  of  problem  on  the  basis  of  the  methods  presented 
in  Reference  1,  that  is  equating  the  stress  level  for  all  failure  modes, 
requires  a  knowledge  of  the  stress  distribution  throughout  the  panel  for 
the  given  loading.  There  are  no  available  stress  analysis  methods,  of 
sufficient  sophistication,  for  corrugated  panels  under  lateral  loading, 
hence  a  large  portion  of  this  report  is  devoted  to  the  development  of  such 
a  method  of  analysis.  The  method  is  essentially  a  deformation  method, 
the  basic  characteristics  of  which  are  described  in  Reference  2j  however, 
a  technique  of  using  Fourier  type  transformations  makes  the  solution  more 
tractable. 

Chapter  1  presents  in  detail  the  method  of  stress  analysis  of  the 
sandwich  panel  with  a  truss  core.  Chapter  2  discusses  the  optimization 
procedure.  In  Appendix  1  it  is  demonstrated  that  the  finite,  one- 
dimensional  structural  element  used  in  the  analysis  is  valid  for  the  type 
of  panel  which  is  considered. 
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NOTATION 


hc 

K 

K*- 

Ai 

n 


Matrix  defined  by  Equation  (1.36) 

Panel  dimension  in  the  y  direction,  inches 
Panel  dimension  in  the  x  direction,  inches 
Stiffness  matrix  defined  by  Equation  (1.2j8) 

Column  matrix  of  Joint  displacements  referred  to 
i-th  coordinate  system  (i  *  1,  2,  3),  inches 
Displacement  of  the  i-th  Joint,  inches 
Number  of  plate  elements  in  panel 
Modulus  of  elasticity,  lbs. /in. 

External  force  per  unit  distance  in  x-direction  applied  to 
i-th  Joint,  lb. /in. 

Internal  force  per  unit  distance  in  x-direction  applied  to 
i-th  Joint,  by  the  J-th  element,  lb, /in. 

2 

Shear  modulus,  lbs, /in. 

Column  matrix  of  Joint  forces  defined  by  Equation  (l.lt3) 

Column  matrix  of  applied  forces  defined  by  Equation  (l*li7) 
in  the  c  ■  2  coordinate  system 

Core  depth,  inches 

(R)(2C)(RT).  See  Equation  (1.55) 

Shape  factor  (dimensionless) 

Width  of  a  plate  element  [i  •  c  (core)  or  f  (face)J  ,  inches 
Number  of  joints  in  panel 

'.'»i.iponent3  of  referred  to  coordinate  system  c 

«here  c  >  1,  ?,  or  ),  lb:;./ir.. 


a 

*s 


Components  of  referred  to  coordinate  system  c,  lbs/in. 
Matrix  defined  by  Equation  (1.50) 

Column  matrix  of  internal  forces  defined  by  Equation  (l.iiO) 
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NOTATION 


'/K  Matrix  defined  by  Equation  (1.36) 

CL  Panel  dimension  in  the  y  direction,  inches 

b  Panel  dimension  in  the  x  direction,  inches 

i  _ 

G  Stiffness  matrix  defined  by  Equation  (1.1*8) 

t 

Column  matrix  of  Joint  displacements  referred  to 
i-th  coordinate  system  (i  *  1,  2,  3),  inches 
Displacement  of  the  i-th  Joint,  inches 
C  Number  of  plate  elements  in  panel 

p 

E.  Modulus  of  elasticity,  lbs. /in. 

External  force  per  unit  distance  in  x-direction  applied  to 
i-th  joint,  lb. /in. 

(  f;  1  Internal  force  per  unit  distance  in  x-direction  applied  to 
i-th  joint,  by  the  j-th  element,  lb, /in. 

2 

G  Shear  modulus,  lbs. /in. 

t 

G*  Column  matrix  of  joint  forces  defined  by  Equation  (l.li3) 

z 

GB  Column  matrix  of  applied  forces  defined  by  Equation  (1.1*7) 

in  the  c  *  2  coordinate  system 

he  Core  depth,  inches 

K  "  (R)(2C)(RT),  See  Equation  (1.55) 

Kr  Shape  factor  (dimensionless) 

X.L  Width  of  a  plate  element  [i  ■  c  (core)  or  f  (face)J  ,  inches 
Number  of  joints  in  pnnrl 

c  P*  QV,*^V  '•''•'•ponent3  of  referred  to  coordinate  system  c 

where  c  a  1,  2,  or  3,  lbr,./ir., 

(  ^  t  i  C  4  r  ^  \ 

pj ,  ,  •’j  Components  of  referred  to  coordinate  system  c,  lbs /in. 

R,  Matrix  defined  by  Equation  ( 1 . 50) 

S  Column  matrix  of  internal  forces  defined  by  Equation  (1.1*0) 
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UAHLR  1 

aVALYSJ?  OF  TRUSS-COPE  SANDWICH  PANEL 


A.  Introduction 

The  analysis  of  a  complex  structure  under  loading  begins  with  the 
decomposition  of  the  structure  into  "basic  elements"  and  assuming  that 
their  response  to  any  loading  that  they  are  required  to  sustain  is 
known.  Viewed  in  this  light,  the  problem  of  determining  the  behavior 
of  the  structure  under  a  specified  loading  condition,  is  really  a 
problem  of  determining  what  loading  will  be  placed  on  each  "basic 
element"  (i,e.  internal  loading)  as  a  result  of  the  specified  loading 
being  placed  on  the  structure  (i.e.  external  loading). 

The  determination  of  distribution  of  internal  loading  for  a  given 
external  loading  condition  can  always  be  accomplished  by  invoking  two 
obvious  physical  requirements.  Tho  first  is  that  every  element  of  the 
structure  will  remain  in  equilibrium.  The  second  is  that  deformations 
of  the  basic  element  as  a  consoquonce  of  the  loading  placed  on  them  will 
not  violate  the  basic  integrity  of  the  structure. 

In  implementing  this  approach  to  the  analysis,  it  is  necessary  to 
choose  the  basic  elements  of  the  structural  system, postulate  the  form 
of  the  loading  that  must  be  sustained  and  determine  the  response  of  the 
elements  to  the  loading.  Great  accuracy  in  this  type  of  finite  element 
analysis  usually  requires  that  the  basic  elements  be  as  3mall  as  possible 
resulting  in  a  structure  consisting  of  a  great  number  of  elements. 
Howover,  the  requirement  that  the  mathematical  formulation  is 
tractable  may  require  that  the  structure  be  divided  into  a  fewer  number 
of  largei  basic  elements. 
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figure  i 

Triangulated  Core  Sandwich  Plate 


It  is  desired  to  determine  the  response  of  this  panel  to  a  lateral 
loading  when  the  panel  is  supported  continuously  along  the  four  edge 


boundaries 
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Due  to  the  construction  of  a  truss-core  sandwich  panel,  it  is  not 
possible  to  easily  determine  the  orthotropic  properties  in  the  x  and  y 
directions  for  purposes  of  analysis  and  for  determining  stress  distribu¬ 
tions  under  lateral  loading.  Therefore,  the  formulation  of  the  problem 
will  be  in  terms  of  more  fundamental  structural  elements  of  the  system. 
One  such  element  is  a  plate  element  shown  in  Figure  ?. 


It  is  assumed  that  for  the  constructions  considered 

i;  «  Xi  «  b  =  0  (X)  i  -  f-  or-  c 

where  A  is  a  characteristic  length  which  gives  a  measure  of  the  rate  of 
spatial  variation  of  all  loading  placed  on  these  elements,  thus  it  is 
possible  to  treat  the  elements  as  one-dimensional.  This  approach,  in 
effect,  reduces  the  problem  to  a  group  of  coupled  one  dimensional 
problems. 

Now  consider  this  formulation  in  greater  detail  when  the  basic 
elements  are  chosen  as  described  above.  A  second  assumption  is  made  that 
t^  is  so  small  relative  to  all  other  pertinent  dimensions  that  the 
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resistance  of  those  elements  to  any  deformation  which  is  out  of  the  plane 
defined  by  the  element  is  neglected.  Consistent  with  this  last 
assumption,  it  is  further  assumed  that  the  basic  plate  elements  are 
loaded  only  by  in-plane  forces  distributed  along  their  edges. 

Finally,  the  applied  lateral  load  on  the  panel  is  replaced  by  statically 
equivalent  line  loads  in  the  x-directlon  at  the  locations  where  the  web 
and  face  elements  intersect. 

C.  Coordinate  Systems  and  Notation 

In  Figure  3  is  shown:  (l)  a  numbering  system  which  is  to  be  used 

to  identify  the  elements  and  Joints;  and  (2)  three  coordinate  systems 

which  are  convenient  to  use  at  different  stages  of  the  analysis.  The 

x-y  (c  "  1,  2,  3)  axes  are  parallel  to  the  principal  moment  of 
c  c 

inertia  axes  of  the  cross-sectional  areas,  also  the  x  axis  is  parallel 
to  the  flute  direction  of  the  core  and  y^  is  parallel  to  the  1 
dimension  of  the  plate  elements. 


Figure  3 

Coordinate  Systems 
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Referring  to  Figure  3*  the  following  notation  is  introduced? 


e  ~  number  of  plate  elements  in  panel 

h  s  number  of  joints  in  panel 

^F*]  =  external  force  per  unit  distance  in  x  direction 
applied  to  the  i-th  joint 


'pI'qVr* 

tFil 


components  of  referred  to  coordinate  system  c 

where  c  r  1 ,  2  or  3 

Internal  force  per  unit  distance  in  x  direction 
applied  to  i-th  joint  by  J-ih  element 

t  i  c  i  c  i  ci. 

Pi .  Si'  ri  components  of  £fj  ^  referred  to  coordinate  system  c 

^•is  K  Cl  rt  1 


"i-T/  * ’Pi 

c ,•  K  c-  t 

ft  *  Pi  -  ‘Pi 
°3  *Yi  *‘3. ‘i 

3/  *V-V 


(1.0  a) 


define  the  internal  forces  acting  on  J-th  element.  Note  that  i  and  k 
refers  to  the  two  Joints  defining  J  with  k  >  J.  The  subscript  J  on  c 
denotes  that  the  coordinate  system  to  which ff^}  is  referred  is  chosen 
so  that  the  y  axis  is  parallel  to  the  direction  of  the  J~th  element. 


W  *  displacement  of  the  i-th  Joint 
cm\VY  ur1--  components  of  referred  to  coordinate  system  c 


oC, i  ‘)i 

py  =Cv  -  V)il 

V;  Pit' Sir  jr 

Sj  iV‘-V‘)t 


(1.0  b) 


define  the  state  of  deformation  of  J-th  element.  Again  i  and  k  refer 
to  the  two  Joints  defining  J  with  k>i.  The  subscript  J  on  c  has 
the  same  meaning  as  above. 
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D.  Formulation  of  Problem 

The  problem  to  be  solved  is  stated  as  follows?  Given  ths 
etructure  and  either  the  external  force  per  unit  distance  to  be  applied 
to  each  Joint  or  the  resulting  displacement  to  be  obtained  at  each  Joint, 
find  the  resulting  internal  force  distribution.  This  report  will  con¬ 
centrate  on  solving  the  problem  when  the  external  forces  are  prescribed 
at  all  Joints. 

It  was  pointed  out  in  the  introduction  that  the  solution  is  obtained 
hy  invoking  the  physical  requirements  of  equilibrium  of  all  points  in  the 
etructure  and  ensuring  continuity  (maintaining  structural  integrity) 
throughout  the  structure  under  loading. 

(1)  Equilibrium 

The  equilibrium  requirement  applied  to  each  of  the  joints 
results  in  the  following  set  of  algebraic  equations  which  expressed  in 
vector  form  are 


ho 


l  =  l|"',h  (It  1) 


In  these  equations  the  summation  is  carried  out  over  all  elements  (j) 
which  frame  into  the  i-th  Joint.  Equations  (1,1)  constitute  a  set  of 
linear  algebraic  equations  on  the  unknowns  .  To  ascertain  whether 

they  are  sufficient  to  determine  the  unknowns"  requires  us  to  simply 
count  equations  and  unknowns.  There  are  obviously  3n  scalar 
equations,  and  since  there  are  e  (•  2n  -  3)  elements,  there  are 
2*  (■  lin  -  6)  vector  unknowns.  Note  that  by  virtue  of  the  assumption 
that  no  elements  can  sustain  a  force  perpendicular  to  the  plane  defined 


See  Reference  2, 
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by  the  element,  each  of  the  vectors|f ^ jia  confined  to  be  in  a  known 
plane, hence  only  two  scalar  unknowns  are  associated  with  each[fj.} 

The  total  number  of  unknowns,  therefore,  is  he  •  8n  -  12  which  is  more 
than  the  number  of  equations  if  n>2. 

Referring  these  n  vector  i^quationa  (1.1)  to  coordinate 
system  c  *  2,  the  following  }n  scalar  equations  results 


+"  2P*  25  0 

i  »  / ,  ■  •  • 

(1.1  ») 

4-V  *0 

L*  1,  •  *  •  n 

(1.1  b) 

Z'rj 

JO 

M 

O 

i  *  /,  *  *  •  h 

(l.l  c) 

(2)  Internal  Force  ~  Deformation  Relationships 

Of  all  of  the  sets  of  unknowns  that  will  satisfy  the 
equilibrium  requirement  expressed  by  Equations  (1.1),  the  one  that 
will  represent  the  actual  distribution  of  internal  forces  in  the 
structure  gives  rise  to  a  deformation  for  each  of  the  basic  elements 
that  will  not  violate  consistent  deformations  of  the  structure.  To 
apply  such  conditions  the  deformation  of  each  element  in  terms  of  the 
internal  forces  acting  on  it  must  be  derived. 

Consider  the  J-th  element  together  with  the  associated 
coordinate  system,  relevant  dimensions  and  th8  loads  to  which  it  is 
subjected  as  shown  in  Figure  lia. 
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Figure  ha 


Due  to  the  assumption  of  a  linear  system,  it  may  be  taken  as  the 
superposition  of  the  following  three  problems  shown  in  Figures 
lib,  he,  and  hd. 


where  t 


Figure  hb 

CTU)  a  -1*  '2/  (*)  j*  I,*  ••,<-an-3(1.2) 

p,U)  - ° 4  p  C x.)  -C'pj  (*)  j*  /„•••.<•  2o-3U*3) 
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where  Cj  fe  C.  c 

7*(k)  ~  M  ~  (*)  y«/  -"<?*2o-3(l.li) 


m. 


*•  x 


t r  <>0 

1  4 


where 


Figure  l»d 

TT(x)  =  C;p/W  +  CJp-M  j  %"*■»#.  1.5) 


Subject  to  the  restriction  that  «.  <£)  and  where 

A  is  some  characteristic  length  associated  with  the  loading,  it  is 
permissable  to  analyze  each  of  these  problems  by  a  suitable  one 
dimensional  theory.  The  Justification  of  such  a  step  together  with 
a  derivation  of  the  appropriate  one  dimensional  theory  is  given  in 
Appendix  I.  A  summary  is  given  below: 

Problem  I  (Fig.  kb)  gives  rise  to  a  deformation  characterized 
by  a  displacement  in  the  y^  direction  which  does  not  vary  with  y  and 
a  displacement  in  the  x  direction  which  varies  linearly  with  y  ^ 

[j..e*  displacement  -  ^3  j  (x)^J  .  The  relationship  between 

(^(V'Oand  YjM ,|3#(a)»  as  derived  in  Appendix  I,  is  given  by: 


9 


-2-K1  G;  t-  i  [p,  <o]  -  k‘G,  £j  tj  £  [  U)  j  o.m 


-  e  'tn-'J 


4K&Ai  ^4  dijf Q  (x)l  +2.  KG-t:  d.jy  (x)|  0. 

X;  3(»  -u)  d  Xjl  '*  j  d  JdxL' 

.  J  r  --  .e  «  2‘ 


7) 


_  y  ■»  /,  ••  ,e«2n-3 

Tn  Fqu-if  i  ons  (1.6)  and  ().7)  G-  and  )/'  refer  to  the 


shear  modulus  and  th<*  Poisson's  ratio  of  tho  j-th  element  whereas 
and  tj  are  the  height,  and  thickness  of  the  J-th  element 
respectively.  K ^represents  a  factor  to  account  for  the  shape  of  the 
cross-section. 

Problem  II  (Fig,  he)  gives  rise  to  a  deformation  which  is 
characterized  by  a  displacement  in  the  y^  direction  which  varies 
linearly  with  yj  j^i.e.  displacement  p  &J  (*) j  •  The  relationship 
between  Jj'(jk)  and  Sj  (x)  ,  as  derived  in  Appendix  I,  is 


Problem  III  (Fig.  hd)  gives  rise  to  a  deformation  which  is 
characterized  by  a  displacement  in  the  x  direction  which  does  not  vary 
with  yj.  The  relationship  between  TV.  GO  and  o<^  (x)  ,  as  derived  in 
Appendix  I,  is 


(1.8) 


j  •  •■,e-2o-3  (1-9) 


(3)  Consistent  Deformations 

In  Section  C,  internal  deformations  are  defined  that 
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are  associated  with  each  distribution  of  internal  forces.  In  order  to 
maintain  the  integrity  of  the  structure  only  those  states  of  deforma¬ 
tions  which  are  associated  with  a  unique  displacement  for  each  Joint 


will  be  allowed. 

This  requirement  will  be 

met  if  OCj  ,^j  ,  Y^’  and 

O*  *  C  *  ft 

Sj  can  be  defined  in  terms  of  JL4.t'  ,  *  U. 

,C/  IT1  and  C'V* 

according  to: 

°<-j  ^ 

?  u*  + Cj  u;H 

j  =  I,  •  •  •  •  ,e»*m-3(i.io) 

=( 

c‘  u*  -  c' 

]*•»**•  ■,€®2n-3(i.xi) 

<c  •  y  c:  ;\  t 

/  t r*  +  '  vt‘)t 

t 

j  -  1,  •  *  **,€=20-3  (l.l?) 

Sj-I 

^  -Cj  ctOt 

j  =1,*  *  ’  *,€  =  2n-3(i.i3) 

In  Equations(l.lO)  through(1.13),  i  and  k  refer  to  the  Joints  at  either 
end  of  the  element  J  and  refers  to  the  coordinate  system  with 
x,  y  plane  parallel  to  plane  defined  by  element.  See  Figure  lia. 

Each  of  the  equations  in  the  system  (1.1)  through  (1.13) 
is  either  a  linear  algebraic  equation  or  a  linear  ordinary  differential 
equation  which  are  obviously  independent  of  each  other.  To  show, 
therefore,  that  the  system  is  sufficient  to  uniquely  determine  all  of 
the  unknowns  again  requires  a  count  of  available  equations  and  unknowns. 
It  has  already  been  shown  that  Equation  (1.1)  represents  3n  equations 
and  6n  -  12  unknowns.  Equations  (1.2)  through  (1.5)  give  8n  -  12  more 
equations  but  also  Introduces  8n  -  12  more  unknowns  £i.e. 
andT Tj  \  j  =  1, ♦  •  •  e»2n-3|.  Equations  (1.6)  through  (1.9)  give  8n  -  12 
more  equations  and  also  8n  -  12  more  unknowns  ^i.e.  oCj  ,  ,  Tj  ,and 
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j*  •,*%€»  2.n  - "b  j  .  Finally,  Equations  (1.10)  through  (1.13)  give 
8n  -  12  additional  equations  while  introducing  3n  unknowns  j\.e.  the  3 
components  of  the  n  displacement  vectors  ^d*^  ,  i  -  1*',  n^j.  Notice 
that  {d*^  (i  r  !,•••,  ri)  completely  determines  the  right  hand  3ides  of 
Equations  (1.10)  through  (1.13).  Adding  equations  and  unknowns  gives 
(27n  -  36)  of  each,  hence, the  system  possesses  a  unique  solution. 

E,  Methods  for  Solving  Equations 

There  are  several  ways  for  solving  the  above  system  of  equations. 
Perhaps  the  most  direct  is  to  use  Equations  (1.2)  through  (1.5)  to 
obtain  ^fj^  in  terms  of  OT* ,  (  ,  and  TT^-  .  Then  by  direct 

substitution  of  Equations  (1.6)  through  (1.9)  into  the  result,  we  obtain 
£fj^  in  terms  of  oC j  ,  Tfj  ,  and  .  Next  we  use  Equations  (1.10) 
through  (1.13)  to  obtain  |f^|  in  terms  of  {  d*^  .  Finally, 
substituting  this  result  into  Equation  (1.1)  will  give  n  vector 
equations  on  the  n  unknown  displacement  vectors  [d1]  (i  ®  l,***,  n). 

Once  this  system  of  equations  has  been  solved  and  the  displacement  of 
each  of  the  Joints  obtained,  then,  ail  of  the  other  unknowns  are 
obtainable  by  direct  calculations  as  dictated  by  Equations  (1,2) 
through  (1.13). 

Without  actually  carrying  out  the  above  substitutions  in  detail 
it  is  readily  apparent  that  the  final  system  of  equations  (3n  scalar 
equations)  constitute  a  system  of  linear  ordinary  differential  equations 
with  constant  coefficients.  The  theory  for  handling  such  systems  of 
equations  is  well  defined  and  the  method  of  attack  is  straight¬ 
forward.  In  the  present  problem,  however,  the  straight-forward 
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procedure  would  become  extremely  complicated  by  virtue  of  t.he  large 
number  of  equations  to  be  handled.  Undoubtedly,  the  large  number 
of  manipulations  that  will  eventually  be  required  will  necessitate 
the  use  of  a  high  spend  digital  computer.  It  is  necessary  ,  there¬ 
fore,  that  the  differential  equations  be  transformed  into  algebraic 
equations. 

The  system  of  n  linear  ordinary  differential  equations  may  be 
transformed  into  an  infinite  system  of  linear  algebraic  equations  in 
the  following  manner.  Py  introducing  a  set  (or  sets)  of  functions  which 
are  complete  over  the  Interval  04  X4/3  and  then  taking  the  unknown 
solutions  as  linear  combinations  of  the  members  of  one  or  another  of  the 
sets,  the  problem  transforms  from  that  of  finding  a  set  of  unknown 
functions  of  x  to  one  of  finding  the  coefficients  iri  the  above 
mentioned  linear  sums.  The  equations  governing  these  coefficients 
may  be  obtained  by  direct  substitution  of  a  genetic  form  of  the  linear  sum. 
The  result  will  be,  in  general,  an  infinite  system  of  linear  algebraic 
oquations.  Thus  we  have  reduced  the  problem  to  solving  algebraic  equations 
but  have  increased  the  number  of  equations  to  an  infinite  number  which 
represents  no  advantage  unless  the  number  of  equations  can  be  truncated 
in  some  way  or  unless  the  equations  can  be  uncoupled  from  on«  another  either 
singly  or  in  finite  blocks.  For  the  present  problem  we  can  introduce 
sets  of  functions  for  which  the  latter  occurs  provided  we  are  willing  to 
let  the  problem  dictate  to  some  extent  the  boundary  condition  at  the 
ends  x  °  o  and  x  *  b.  Since,  in  practice,  the  boundary  conditions 
are  not  clearly  defined,  for  example,  clamped  or  hinged,  this  last 
restriction  is  not  very  severe.  See  further  discussion  in  Appendix  I. 
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F,  Transformation  of  Equations 

The  first  step  is  to  Introduce  the  following  sets  of  functions 
which  are  complete  and  orthogonal  over  the  interval  O  ^  X  i  b 

Ism  COm  Y.\ 


and 


^  cos  CO m 


where  CO»n  «  !3Q_2f 
b 


m  «-  \t'Z,  •  •  ,  CO 


Next  we  expand  all  of  the  variables  in  Equations  (1.1)  through 
(1.13)'  in  terms  of  one  or  the  other  of  the  above  sets  according  to 
Table  1. 

Table  1 

Cosine 


Sine 

« 

y  and  z  components  of 
y  and  z  components  of  ^  fj  j 

*; 

y  and  *  components  of 

Introducing  the  notation 
00 


x  component  of 
x  component  of  £  f  j} 

Pi 

Pi 

cc;  ^ 

x  component  of 


Sin  COm  X 

OP 

COS  OOrnX 


By  substitution  of  Equations  ( 1 . li» )  and  (l.ljj)  into  Equations  (1.1) 
through  (1.13)  and  making  use  of  orthogonality  when  appropriate,  the 
following  system  of  algebraic  equations  are  obtained.  Equation  (1.1) 

1U 


(1.1b) 

(1.15) 


r.w:-ASL-lUl 

is  referred  to  coordinate  sysVm  c  -  2,  which  is  taken  na  the  global 
reference  exes. 


z^jt  +  cp*-):  -  o 

(1.16) 

00  «  q,  •  •  ft  co 

j 

| 

1 

Z('a;;)MlQ;):-o 

n 

(1.17) 

W  "  O,  •  •  ft  00 

+ 

1 

j 

1 

Z(v;)>m;):--o 

t  «■  /  •  • 

(1.18) 

-  O,  •  •  *,oo 

! 

(CTj)Z,  =  C'^-)l+  {%}  )l 

j  *  <r  *  ,  e 

(1.19) 

•on  *  O,  *  '  ftOO 

l 

j 

(p,)°  -  C'p-)L-{c,p‘Ym 

J  *•>•••>  e 

(l.?0) 

t'V'  ’  O,  •  •  ftOO 

1 

! 

j  -  i,*  •,  e 

^  (i-?1) 

nr»  -*0,  ••  ft  co 

i 

i 

1 

i 

(rr);-C'p;)>r'p;): 

j  -  l*  •  ••»  e 

(1.??) 

^  -  o,  •  •  }00 

l 

j! 

1 

j 

1 

i 

] 

(cr):  - [2 K >6,1,  w„]  (ft);  .k' G,  ft  ft 

( l.?3) 

J 

I 

i 

] 

] 

l 

J-o  -e 

^=0,  •  ',  CD 

.  «  1  .  .  .  £» 

I 

i 

i 

1 

t>  ' 

0  -O,.  •  ■  ,co 

J  ’’  *!  (l.?$) 

rv  *  0#  •*,  CO 

j 

| 

k>: 

/  '  . 

J  =  *'  "  e  (1.26) 
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| 
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[(  'u.*)~,  +  (  '  U.‘ )„ 


((5j)l-z[C'u‘)l -Cuc')l] 


(^•)H|V)>(V), 


s 

rv-\  I 


**,e 

nr>  -  O, 

*  •*,  00 

* 

u- 

nr\  - 

0, 

***,00 

»>• 

•  *>  e 

rr\  - 

=  o, 

•  •  *,oo 

J  s 

V 

*•>  e 

rr»  - 

0, 

•  • 00 

Equations  (1,16)  to  (1.30)  are  the  equations  governing  the 
coefficients  in  the  expansions  of  our  unknown  solution  functions.  As 
before  in  dealing  with  equations  (1.1)  to  (1.13)  the  above  equations 
can  be  solved  by  taking  as  the  basic  unknowns  the  components  of  the 
displacements  of  the  Joints (i  "  1,'",  n),  solving  for  them,  and 
then  obtaining  all  other  unknowns  by  direct  calculation.  Referred  to 
the  c  *  2  coordinate  system,  the  basic  unknowns  are  the  functions 

, 1  \fL  (X  )  f  and  1  US' 1  (x. )  ^4  =  /,  nj  or  the  coefficients 
CU1)"  f  and  £*//•  c^Tho 

procedure  is:  (1)  substitute  the  latter  coefficients  into  Equations 
(1,27)  to  (1.30)  and  introducing  a  coordinate  transformation  when 
necessary)  (2)  substitute  the  results  into  Equations  (1.23)  to  (1.26)) 
(3)  substitute  these  results  into  Equations  (1.19)  to  (1.22))  and 
(li)  finally  substitute  these  results  into  Equations  (1.16)  to  (1.10) 


(l.?7) 


(1.28) 


(l.?9) 


(1.30) 
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to  obtain  an  inf'nKe  sy-t**m  of  algebraic  equations  on  the  infinite 
number  of  unknown:'.  Koto,  however,  that  the  equations  governing 
OU4)*^  ,  (  *"  (/*')«*  ,  «nd  for  one  particular 

value  of  m  are  not  coupled  to  those  for  a  different  value  of  m. 

This  means  that  the  infinite  number  of  equations  uncouple  in  blocks 
of  3n  linear  algebraic  equations  on  the  3n  unknowns  ^i.e,  (lU.*’)rv>> 

(  f  ),y,  ,  C  {~K-f  )e-\  [_t  —  /,  •  «  -  n  *  rr\  -  SO*V~\€  irVt"«gC*~J^ 

and  the  problem  reduces  to  a  series  of  problems  involving  3n  linear 
algebraic  equations.  Of  course,  for  an  exact  answer,  we  must  solve 
an  infinite  number  of  such  problems.  For  practical  purposes,  however, 
only  the  first  few  coefficients  will  offer  a  significant  contribution, 
hence,  solving  only  a  few  problems  will  give  a  sufficiently  accurate 
answer. 

The  question  of  the  boundary  conditions  at  x  »  0  and  x  -  b  has 
to  be  answered.  Tho  sets  of  functions  used  in  our  expansions  will  result 

in  y-U).  iei’>  =  4^±>=  O 

*  d*  d  x. 

regardless  of  the  values  of  the  coefficients  obtained  in  solving  the 
above  set  of  equations.  This  condition  corresponds  to  simple  supports 
for  the  flexure  deformation  of  tho  basic  elements;  freo  ends  for 
extension  in  the  x  direction  for  tho  basic  elements  and  no  extension 
in  the  y  (width)  direction  for  the  basic  elements.  Although  the  actual 
conditions  may  be  slightly  different,  the  error  introduced  by  assuming 
these  conditions  should  bo  confined  to  the  immediate  vicinity  of  the 
supports.  Notice  that  other  boundary  conditions  could  be  treated  but 
that  would  require  placing  certain  restrictions  on  the  values  of  the 
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coefficients  that  are  to  be  obtained.  These  restrictions  would  destroy 
the  uncoupling  achieved  above  which  would  necessitate  some  type  of 
iteration  scheme  for  &  solution. 


0.  Matrix  Formulation  for  m-th  Coefficient 

A  matrix  formulation  is  presented  now  to  obtain  the  solution  of  the 
system  of  equations  involving  the  coefficients  of  the  infinite  series 
introduced  in  Section  F. 

As  indicated  in  Section  F,  the  basic  unknowns  are  the  appropriate 
transform  of  the  components  of  {d*^  (i  *  n)  referred  to  the  c  ■  2 

coordinate  system.  Arranged  in  a  column  matrix  of  }n  elements,  the 
unknowns  to  be  determined  are 


x 


D 


c*to: 

cu'): 

e 

C  'ir% 
(V% 
« 

(VI 

(v): 

(V£ 

• 

(Vi 


a. 31) 


The  first  step  is  to  introduce 
system  to  the  c  ■  l  system  and  the  c 


transformation  matrices  fro*  the  c 
■  3  system. 


•  2 
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and 


whore : 


~[l] 

[o] 

[°] 

[0] 

cos  e[l] 

-  sir\  e[l] 

[o] 

S  in  ©[  ]  J 

COS  ^[i] 

[I] 

[°] 

[o] 

[o] 

cos9[l] 

Sin©[l] 

[0] 

-5in  e[j] 

CO*  ©  [i  J 

_ 

(1.32) 


O  is  shown  on  Figure  1 
to]  is  the  n  x  n  null  matrix 
[i]  is  the  n  x  n  unit  matrix. 

The  appropriate  transformation  equations,  as  becomes  apparent 
by  expansion,  are 


1 

D 


T.t  D 


(1.33) 


D  -  Tai  D 

and  are  defined  by  the  matrix  given  in  Equation  (1.31)  except  the 
scalar  elements  are  given  with  resnerh  to  the  c  *  1  and  c  ■  3  axes 
respectively. 

To  represent  Equations  (1.??)  through  (1,30)  in  matrix  form,  sub¬ 
divide  the  elements  of  the  structure  into  three  groups  according  to 
which  coordinate  system  has  the  y  -  x  plane  parallel  to  the  plane  of  the 
element.  Referring  to  Figure  3,  the  groups  have  the  following  numbers. 
Oroup  li  1,  5*  9,  *  *  •  ,  e  -  2 
Group  2t  2,  It,  6,  ...,e-l 
Group  3»  3,  7,  *  *  •  *  ,  e 
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e  »  1  „  n  -  1 

Groups  1  and  3  each  have  ^  ?  members  whereas  Croup  2  has 

e  “  1 

2  *  n  -  2  members.  Now  for  each  of  the  above  groups  a 

rectangular  matrix  is  constructed  with  the  number  of  rows  equal  to 
the  number  of  members  in  the  group  and  the  number  of  columns  equal  to  the 
number  of  joints  in  the  structure  according  to  the  following  rule. 
Associate  each  row  with  a  member  of  the  group  (i.e.  first  row  with 
element  1,  second  row  with  element  etc.)  and  each  column  with  a 
joint  of  the  structure.  Now  we  set  every  matrix  element  in  the  row 
associated  with  a  given  structural  element  equal  to  zero  except  the  two 
which  correspond  to  the  joints  into  which  the  structural  element  is 
framed.  These  are  set  equal  to  +1  or  -1  according  to  whether  the  Joint 
number  is  larger  or  smaller  respectively.  Restricting  our  attention  to 
Group  1,  the  matrix,  so  defined,  is 


-1  1  0  0  0  ....  0 

0  0-1  1 . 0 

0  0  0  0  -1  ...  .  0 

'  a  ■  • 

»  «  « 

9  •  f 

0  0 . 1  1  0 

With  the  aid  of  this  matrix  ths  equations  (1.27)  through 
(1.30)  for  elements  that  fall  within  the  first  group  may  be 
written  as 

‘E  =[VD]i; 


(1.310 
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r  I/- . )  - 

C«.s) 


vrhere 


C. 

nr\ 


E  = 


•  *  N*- 

(.(Xe-iAv 

cp»): 


*  \e 
(  (3  «.-l  jm 

cy.x: 

c  r.^ 

c  i<-X. 
<  b,U 
(Ss)^ 


nr\ 

a 


U.A 


and 


'A  = 


’a 


[°: 

[° 


[o]  [*a]  [o] 

L[0]  M  [qj  J 


(1.35) 


(1.36) 


in  which  [d]  Is  obtained  from  'a.  by  replacing  all  -1  with  a  *1  and 

[o]  1®  an  *fi)null  matrix.  That  matrix  equation  (13b  )  does 

indeed  represent  Equations  (1.27)  through  (1.30)  for  the  numbers  of 

Group  1  may  be  seen  by  direct  expansion.  Similar  results  may  be 

obtained  for  the  other  two  groups  of  members. 

The  deformation  of  the  elements  in  the  three  groups  are  express- 
2 

ible  in  terms  of  the  D  column  matrix  by  Equations  (1.33)  and  (1.3b) 
as  follows i 


lE=‘ATuD  Me  T„  »  [II) 


(1.37) 


To  represent  Equations  (1.23)  through  (1.26)  in  matrix  form  the 
subdivision  introduced  in  the  proceeding  paragraphs  is  kept  but  introduced 
here  is  the  restriction  that  the  properties  of  all  elements  in  Oroups  1  and  3 
are  the  same  and  properties  of  all  elements  in  Group  2  are  the  same.  Under 
this  restriction  the  following  matrices  are  Introduced. 
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where  the  index  i  (i  «  1,  2,  3)  indicates  the  group  of  members  being 
considered  to]  and  Cl  3  are  squire  null  and  unit  matrices  respectively 
with  number  of  rows  and  columns  equal  to  the  number  of  members  of  the 
group.  Consequently,  for  those  elements  of  Croup  i,  the  Equations 
(1.23)  through  (1.26)  may  be  written  as 


'5  =;£;E 


(1.39) 


where 


(rrX 

(rrX 


(1.1*0) 


Note  t  CJ.o^p  i  ,  r=e- Z 
G  ^oup  2. ,  r  -  € 
Group  3,  1“  =  0  —  I 


(aX 
( x, ): 

cm: 

(XX 

Combining  Equations  (1.37)  and  (1.39),  the  internal  forces  in  the 

2 

elements  of  ail  three  groups  ra,iy  be  expressed  in  terms  of  D  according  to 


'S  =  2  ‘A  T, t *  L> 


(  «•  -  I,  2,3,  and  la  =  [i])  (1.1,1) 
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The  final  step  is  to  obtain  the  forces  that  arise  at  each  joint 
doe  to  the  internal  forces  in  the  elements,  sum  all  of  these  forces  and 
substitute  in  Equations  (1.16),  (1.17),  and  (1.18),  This  must  be  done 
separately  for  each  of  the  three  groups  previously  defined  since  each 
group  will  give  the  joint  forces  with  reference  to  a  different  coordinate 
system.  Consider  initially  the  first  group  of  members  then  it  may  be 
shown  by  direct  expansion  that 

G,  =  At  ' 5  U.b2) 

where  *AT  is  the  transpose  of  the  matrix  A  and  G,,  is  a  column  matrix 
representing  the  sum  of  all  forces  acting  on  the  joints  by  elements  in 
the  first  group,  referred  to  the  coordinate  system  c  *  1. 


(1.143) 
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A  similar  expression  may  be  written  for  the  other  two  groups, 
so  that  in  general 

_  i  A  T  i 


Go  =  1 A  5 


(l.Ui) 


where  i  *  coordinate  system 

k  ■  group  number  of  members. 

To  add  the  effrcts  of  the  three  groups  requires  that  all  forces 
be  referred  to  the  same  coordinate  system  tx.e.  c  a  2  system).  This  is 
done  by  using  the  transformation  matrices  anc*  T32  gi^en  in 
Equation  (1.3?).  Por  the  contribution  of  group  k,  we  have 
2  _  , \  T  A 


Gk-CTJ  gk 

where  (TkJ1  is  tho  transpose  of  Ttt  and 


(l.W) 


Adding  together  *0j_,  G^#  and  ^0 ■  ,  the  total  force  applied  at  each 

Joint  by  the  basic  elements  is  obtained  and  Equations  (1.16)  through 
(1,18)  reduce  to 


ZG ,  +  Gi  1  Gs  +  Qc  =  O 


(1.1*6) 


2r~ 

where  Oa  is  the  external  force  applied  referred  to  the  c  -  2  coordinate 


system.  That  is, 


cp): 

(lPX 

(ir): 

cqx 

cq‘>: 

cq”)1 

crx 
c  rx 

cra: 


(1.17) 
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Combining  all  of  the  previous  results  0^,  Gj  and  0^  are 

2 

expressed  in  terms  of  the  single  unknown  matrix  D  and  the  result 
substituted  in  Equation  (I.I16)  to  give  the  matrix  aquation 

*C'D  +  e  -  0  (1.U7) 

where 

lc  -  )(.'f\)iT,i)A'!\)rcz)c/\)  «•  d.1,8) 

(T,x)TA)tCZM’A)(T«) 

is  the  stiffness  matrix.  The  matrix  of  unknown  deflections  can  now  be 

2 

obtained  from  Equation  (l.ii7)  in  terna  of  the  inverse  of  C  and  is 
expressed  as 

2D=-(*C)  (1.1,9) 

Since  all  of  the  desired  unknown  quantities  are  obtainable  from 

2 

D  by  direct  calculation,  Equation  (1.2,9)  represents  the  formal  solution 
to  the  problem. 

Before  proceeding  any  further,  the  matrix  formulation  is  summarized 
in  Figure  It.  The  lines  represent  matrix  products  with  the  multipliers 
being  indicated  over  the  line.  Joining  of  two  or  more  lines  indicates  a 
summation. 
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Figure  U 

Summary  of  Matrix  Computations 

H.  Method  for  Solving  Large  System  of  linear  Equations 

As  already  stated,  the  formal  solution  to  the  problem  is  containod 

in  the  matrix  equation  ( 1 . l»9 ) .  Implicit  in  this  simple  matrix  equation, 

2 

however,  is  an  inversion  of  the  large  order  matrix  C.  Indeed,  in 
carrying  out  the  solution  of  a  given  problem  it  is  numerically  easier 
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not  to  make  use  of  the  formal  solution  of  the  system  of  equations 

represented  by  Equation  (1J*7)  (i.e.  not  t.o  obtain  (*C)~  )  but  to 

obtain  the  solution  in  a  different  nanncr.  In  this  section  an 

algorithm  to  accomplish  this  goal  is  developed. 

The  following  interpretation  can  be  given  an  element  of  the 
2  2 

C  matrix.  C  represents  a  component  of  the  net  internal  force  which 

y 

acts  on  a  joint  which  is  defined  by  i  for  a  displacement  field  which 
consists  of  a  unit  displacement  in  a  direction  and  for  a  joint  defined 
by  while  .1 1  '  ‘it ui  i  nutjons  are  kept  zero.  It  is  obvious,  there¬ 
fore,  that  unless  the  joint  defined  by  j  is  directly  connected  to  the 
Joint  defined  by  i  by  an  element  of  the  structure,  then  must  be 

y 

zero.  This  loads  to  t.bn  ,  .,n  iu*non  that  the  majority  of  elements  of 
2 

the  "Cy.  matrix  are  zero.  The  solution  of  the  equations  represented  by 
( 1.U7 )  can  be  most  easily  accomplished  by  partitioning  into  smaller 
groups  of  coupled  equations. 

Before  accomplishing  this  partitioning  it  would  be  helpful  to  first 

recast  the  formulation  in  such  a  manner  that  all  of  the  non-zero  elements 
2 

of  the  C  matrix  cluster  around  the  main  diagonal.  This  is  done  by 
introducing  the  following  square  matrix  which  contains  only  ones  and 
zeroes. 
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where  the  Kronecker  delta 


I  if 

o  I  f-  l  ?  i 


has  been  introduced.  It  may  be  noticed  that 


RRT  •  I 


f 

where  R  indicates  the  transpose  of  R  and  I  indicates  the  3n  x  3n 
unit  matrix. 

Multiplying  Equation  (1.1*7)  by  R  and  noting  Equation  (1.52) 
results  in 


(R2CRT)  (R)(2D)  -  -  (R)(20g) 


(l.5o) 


(1.5D 

(1.52) 

(1.53) 
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where, 


(rX’d)= 


CP')~ 

CtrX 

CQX 

CR')1 

cu‘): 

CPl)t 

(v*>: 

« 

=X  (r)Cg=)= 

CQ')l 

(’  R’)m 

9 

cu" ): 

(*  PX 

Cr'K 

CQ'K 

C  R")£ 

=Y 


the  new  stiffness  matrix 

K  -  (R)(2C)(R  f 


(1.5M 


U.55) 


now  has  all  non-zero  elements  clustered  around  the  main  diagonal.  The 
system  of  equations  expressed  in  the  form  given  in  Equation  (1.53)  is 
solved  rather  than  the  fo**r>  f»ivon  in  Equation  (i.li7). 

Expanding  the  matrix  Equation  (1.53),  the  resulting  system  of 


equations  may  be  represented  by  the  following  partitioned  formt 

Kii  X i  ♦  KitX»  +• 

Ki,  y,  4-  K*sX» 

K  »x  Xa  4-  KhXj  4-  Ks4 


In  Equation  (1,56),  K^, 


KitnlcXe  +  Kee  Xte*0  4-  K«  («*i)  XGnz)  = 
K  U  *1 X  t  I )  x«e  4- 1 )  +  K  c  e  ♦  zXc + 1)  X  (£+*) =» 


represent  submatrices  obtained  by 


partitioning  the  K,  X,  Y  matrices  according  to 


y» 

Yx 

y3  U.56) 
y<fe*-o 


K 


Kh  Kix  *  #  »  ,  , 

'  *  . 

Kti  Kii  Xxi  •  ,  t  •  ,  , 

*  Kji  K*a  Ka^ 

«  »  •  !  Kcp+')e  Kcpv»Xp*0  KlpuXptz) 

Kcp^Xp*-*)  Ktp^Xp+2) 

•  I  •  t  #  # 

(1.57) 
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X  = 


X, 

Xx 

n 

• 

• 

Y= 

,  .  ,X  -C 
»*  - 

_ 1 

# 

Xp*t 

Yprx 

In  Equations  (1.56)  and  (1.57),  p  is  given  by 

/>  -  f  ~  I 


(1.58) 


and  the  submatrices  have  dimensions  as  indicated  in  Table  2. 


Table  2 


Matrix 

Ku 

Ku  (i  -  2,3 ,..(>*1) 

K  (p  +2 )  (  P  +2 ) 

K12 

K2l 

Y.  (i  -  2 ,3,..,p) 
til*  t) 

K  (i  +  l)i  (i  *  2,3, ...,p) 

K  (f>+  l)(p*  2) 

K  (p+  2 )  ( p  1) 

X1 

X1  (i  -  2,  •  •  •  ,  p  ♦  1) 

V  +  2 


Ti 


2»  *** »  P  +  1) 


p+  2 


Dimension 
6x6 
12  x  12 
6x6 

6  x  12 
12  x  6 
12  x  12 
12  x  12 
12  x  6 
6  x  12 
1x6 
1  x  12 
1x6 
1x6 
1  x  12 

1x6 
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The  system  of  equations  expressed  in  the  form  indicated  by  (1,56) 
say  b#  solved  as  follows:  From  the  fH3t  wr  obtain 

X1  •  KU  [li  -  K12  x2 

Using  this  in  the  second  we  can  obtain 

^  *  (k22  -  Ka  K-J  Kurl  (r2  -  x21  x-J  rx  -  k2j  x3). 

This  Is  then  used  in  the  third  to  obtain  an  expression  for  in  terms 

of  X,  which  is  then  used  In  the  fourth  to  obtain  X,  in  terms  of  X  .  etc. 

U  u  5 

f'roceeding  in  this  manner  we  arrive  at  the  (p  *  2)  equation  expressed 

only  in  terms  of  ,  which  may  be  solved.  Once  X  _  has  been  so 
^  P-t 2  P*  2 

determined  X  ,  then  X  then  X  ,  .  ,  .  X  may  all  be  obtained  by 

p*  i  p  r  -i  i 

back  substitution.  Thus,  the  solution  may  be  expressed  by  the 
following  set  of  equations: 

vC(pvi)  -  ©  c  p  ♦  i  )(p  *  i)  2.  t  p  +  z> 

^X-Cpn)  -  Bfp» [Z  tp+-i) -K<p*  'Xp»*)Xcp*i^1*'^ 
X- p  £  Bpp  zLp  ~~  KpCp+-OX.p^j 
X,  0  . ,  j/E.,  -KaXj 


where 


Bn  »  K,,'1 

8*1  =  -  Hvzi  B„ 

Bj3  --  (Kj3-  Ksz  Bn  A«)' 


(1.60) 


B(pnXpn)  (K(pf  jxpvt)  Kto  ez^Xpei)  S(p*iXprt)  K.  (pn  Xp*  z)) 
and 

Z,  =  Y, 

Z,  =  Ya  -  Kn  BuZ,  (1.61) 

ZLa  =  NY  -  K32.B2iZ«i 

ZL<p*0  =  Yp  n)  “  Ktp*a)(p+0  Btp»»'sCp+-i)<Zp4.o 

It  is  important  to  note  that  the  solution  as  expressed  by 
Equations  (1.59)  through  (l.6l)  only  requires  the  inversion  of  a  12  x  12 
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matrix  and  is  independent  of  < he  number  of  joints  in  the  structure  to 
be  analyzed.  It  ought  also  be  not'-d  that  the  programming  of  the 
solution  as  represented  here  is  a  relatively  simple  task  and  that  a 
program  can  be  easily  written  for  a  general  value  for  n. 


I.  Synthesis  of  the  Final  Solution 

Based  on  the  explanations  given  in  Sections  F  and  0,  the  synthesis 

of  the  final  solution  for  the  stress  distribution  in  the  panel  is  now 

summarised.  The  set  of  3n  equations,  essentially  given  by  Equation 

(1.1*9),  are  solved  a  finite  number  of  times  for  the  coefficients  of 

the  transformed  displacement  components,  denoted  in  Equation  (1.31)  by 
2 

D.  This  finite  number  of  solutions  of  Equation  (1.1*9)  will  depend  on 
the  characteristics  of  the  lateral  loading.  However,  for  the  types  of 
flight  loadings  usually  encountered, the  number  of  terms  in  the  aeries 
(or  coefficients  of  the  transform  of  displacements)  represented  by 
Equations  (l.ll*)  and  (1.15)  will  be  probably  less  than  five.  Once  the 
displacement  coefficients  are  determined  then  the  internal  forces  can 
be  obtained  by  using  Equation  (1.1*1),  Finally  the  actual  joint  dis¬ 
placements  and  stress  distribution  in  the  panel  ore  obtained  by  summing 
the  finite  number  of  terms  in  the  series  of  the  form  given  in 
Equations  (l.ll*)  and  (1.15)  and  in  accordance  with  the  scheme  of 
Table  1. 

J.  Remarks 

The  method  of  analysis  developed  in  this  chapter  has  been  based 
on  the  assumption  of  linear  elasticity  and  isotropy.  However,  it  is 
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possible  to  extend  the  technique  t  include  orthotropic  materials  by 
deriving  the  force-displacement  relationships  similar  to  those  given 
in  Equations  (1*6)  through  (1.9)  and  discussed  in  detail  in  Appendix  1. 

The  development  of  the  analysis  algorithm  in  this  chapter  was 
based  on  a  truss-core  panel.  However,  the  results  are  presented  in  the 
fora  of  matrix  formulas  which  are  completely  general.  Hence  any  change 
in  core  geometry  (®,g.,  web  core)  is  reflected  only  in  those  matrices 
which  define  the  geometry  of  articulation,  such  as  Equations  (1.32), 
(1.33),  and  (1.36). 

Although,  beyond  the  scope  of  this  work,  it  appears  advisable  to 
investigate  the  stress  distributions  in  corrugated  core  panels  under 
lateral  loads  by  utilizing  ct  computer  program. 

It  is  recommended  that  several  simple  lateral  distribu¬ 
tions,  for  example  uniform  or  triangular,  be  considered  to  obtain  the 
stress  distributions  Tf  ,  ,  <T  ,  in  various  elements  of  the  panel. 

By  systematically  investigating  the  stress  distributions  as  various 
parameters  are  varied,  such  as  thicknesses,  aspect  ratios,  different 
materials  for  faces  and  core,  it  appears  feasible  that  empirical 
expressions  could  be  formulated  so  that  it  would  not  be  necessary  to 
resort  to  the  more  time  consuming  matrix  analysis.  Further,  once  such 
closed  form  expressions  for  stresses  are  available  the  method  of 
structural  optimization  used  in  Reference  1  can  be  applied. 
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OPTIMUM  DESIGN  OF  TRUSS-CORE  SANDWICH  PANELS 
U*  Introduction 

7*19  first  stop  in  the  development  of  any  optimum  design  procedure 
is  to  precisely  define  what  is  to  be  meant  by  optimum  and  then  to  trans¬ 
late  this  definition  into  a  mathematical  language.  There  are  several 
different  ways  in  which  this  can  be  accomplished  and  two  of  these  are 
discussed  here. 

One  common  method  of  defining  optimization  is  first  to  define  some 
function  or  functional  of  the  design  parameters  and  to  equate  optimiza¬ 
tion  of  design  with  either  the  maximization  or  the  minimization  of  this 
function  or  functional.  Two  frequently  used  choices  for  the  function  to 
be  extremized  are  the  total  cost  of  the  structure  and  the  total  weight 
of  the  structure.  It  is  important  to  emphasize  that  the  design  which  is 
optimum  from  a  total  cost  point  of  view  is  not  necessarily  optimum  from 
a  total  weight  point  of  view  and  vlce-verea.  It  would  only  be  known,  a 
priori  ,  that  both  criteria  would  lead  to  the  same  design  if  it  is  known 
that  the  only  factor  that  went  into  the  total  cost  is  the  weight.  Never¬ 
theless,  the  tendency  quite  often  is  to  choose  the  minimum  weight 
criterion  when  the  desire  is  actually  to  minimize  cost  since  it  is  usually 
much  easier  to  express  the  functional  dependence  of  weight  on  the  design 
parameters  than  it  is  to  express  the  functional  dependence  of  cost  on  the 
design  parameters.  In  this  work,  the  weight  factor  alone  is  considered 
and  all  other  factors  are  ignored,  A  modified  version  of  such  an  approach 
would  be  to  investigate  the  minimum  weight  design  and  determine  if  there 
are  ary  excessive  costs  associated  with  it  (i.e.  possibly  high  fabrication 

3h 


MEC-ASL-1111 


coats).  If  none  appear,  then  minimum  weight  is  classified  as  optimum,  but 
if  excessive  costs  are  uncovered  then  the  structure  i-  modified  to  remove 
the  excessive  cost  problem  and  a  somewhat  higher  weight  structure  is 
classified  as  minimum. 

An  optimization  procedure  as  out  Lined  above  leads  to  a  mathematical 
problem  which  falls  within  the  class  of  problems  termed  extremum 
problems.  The  approaches  for  obtaining  their  solutions  fall  Into  two 
distinct  classifications,  one  is  termed  the  direct  approach  and  the  other 
is  termed  the  indirect  approach.  The  direct  approach  consists  of  simply 
evaluating  the  function  to  be  extremized  for  all  possible  combinations  of 
the  design  parameters  and  then  comparing  the  results  to  see  which  is 
the  desired  extremum.  It  is  apparent  that  such  an  approach  can  easily 
become  intractable  as  the  number  of  design  parameters  Increase  or  the 
spectrum  width  of  values  is  broadened.  On  the  other  hand,  it  is  often 
possible  to  generate  a  scheme  to  "zero  In"  on  the  optimum  value  even  for 
large  numbers  of  design  parameters. 

The  indirect  approach  is  to  develop  some  easily  investigated  require¬ 
ments  that  a  particular  design  must  satisfy  if  it  is  to  extremize  the 
desired  function  and  to  throw  out  all  designs  that  do  not  live  up  to 
these  requirements.  The  designs  that  sat isl'y  these  requirements  are  then 
subjected  to  the  direct  approach.  As  an  example,  it  is  possible  to  show, 
subject  to  certain  continuity  requirements,  that  if  the  function  of  the 
design  parameters  13  to  take  on  an  extreme  value  for  a  particular  design 
then  any  small  change  In  the  design  will  produce  no  change  in  the  function. 
This  fact  can  be  translated  into  mathematical  language  and  the  result 
used  as  a  test.  Although  the  indirect  approach  is  frequently  much  faster 
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than  the  direct  approach  it  i?  not  without  its  problems  since  "easily 
investigated  requirements"  are  frequently  not  easily  investigated  and 
sometimes  no  design  will  pass  the  postulated  requirements  which,  of 
course,  does  not  mean  there  is  no  optimum  design  but  rather  that  the 
requirements  are  too  severe. 

As  a  rule  the  indirect  approach  will  work  quite  well  if  the  design 
parameters  are  not  limited  m  any  way  or  if  the  restrictions  to  be  placed 
on  the  design  parameters  are  expressible  by  equations.  Difficulties 
arise  when  the  restrictions  are  presented  in  expressions  containing 
inequalities.  For  the  design  of  a  structure,  these  latter  restrictions 
usually  occur  since  the  constraints  that  are  to  be  placed  on  the  para¬ 
meters  are  simply  that  the  structure  will  not  fail  in  any  manner  (i.e. 
any  structure  is  acceptable  that  does  not  fail). 

A  second  approach  to  the  entire  problem  of  optimum  design  of 
structures  is  based  on  the  physical  definition  of  an  optimally  designed 
structure  as  one  which  is  made  to  use  all  of  its  load  carrying  ability 
to  the  greatest  extent  possible.  Starting  with  this  definition  of 
optimum,  the  procedure  is  to  enumerate  the  various  possible  modes  of 
failure  and  then  choose  the  design  parameters  such  that  as  many  modes 
as  possible  are  brought  to  incipient  failure.  If  the  condition  of  failure 
for  each  of  the  modes  results  in  a  system  of  linear  algebraic 
equations  on  the  design  parameters  and  if  all  of  the  algebraic  equations 
bo  obtained  are  linearly  independent  then  a  solution  is  assured  since  as 
many  modes,  as  there  are  design  parameters,  can  be  brought  to  incipient 
failure.  If  the  condition  of  failure  for  each  of  the  modes  results  in  a 
non-linear  equation  then  we  do  not  have  such  an  assurance  and  we  must 
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test  to  see  if  it  is  possible.  See  'inference  1.  Several  recent 
investigations  in  non-linear  structural  problems  O,  h)  indicate 
that  solutions  do  converge  to  disuoct  and  real  -alues  of  the  parameters 
involved  and  thus  unique  solutions  do  -p««uH. 

It  i3  possible  to  compare  this  second  approach  of  optimization 
to  the  first  approach  in  the  following  manner. 

Assume  that  minimum  weight  Is  the  optimization  criterion  in  the 
first  approach.  Mathematically, then,  the  problem  is  to  choose  the  design 
parameters  such  that  the  weight  function  is  a  minimum  when  the  design 
parameters  are  so  chosen  that  the  structure  does  not  fail.  As 
previously  stated,  the  restriction  that  the  structure  should  not  fail  is 
expressed  mathematically  as  a  series  of  inequalities,  one  arising  from 
each  of  the  various  possible  failure  modes.  Physical  reasoning  is  now 
introduced  and  it  is  postulated  that  the  minimum  weight  structure  is 
probably  one  for  which  as  many  as  possible  of  these  inequalities  is  just 
satisfied.  The  object,  therefore,  is  to  choose  the  inequalities  that  are 
closest  to  being  violated  (t.e.  as  many  as  possible)  and  design  so  that 
all  are  ju3t  at  the  point  of  being  violated.  It  is,  of  course,  necessary 
to  check  that  all  of  the  remaining  mequal x ties  are  still  satisfied. 

Although  the  latter  approach  is  based  more  on  physical  reasoning 
than  is  the  former,  it  is  undoubtedly  the  easier  approach  to  formulate 
provided  the  various  conditions  of  incipient  failure  are  expressible  in 
closed  form.  Thus,  it  is  a  superior  approach  from  an  engineering  point 
of  view.  For  cases  in  which  the  conditions  of  incipient  failure  are  not 
expressible  in  closed  form  it  appears  that  the  second  approach  can  only 
be  solved  by  a  trial  and  error  procedure.  The  difficulties,  therefore, 
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appear  to  be  the  same  as  those  which  would  be  encountered  in  minimizing 
a  function  by  the  direct  method,  in  fact,  it  appears  that  there  would 
be  sore  difficulties  since  it  is  necessary  t.c  satisfy  more  conditions 
by  trial  ard  error. 

B.  Design  Parameters 

Consider  the  truss-core  sandwich  plate  shown  in  Figure  5.  It  is 
desired  to  optimize  the  design  of  this  structure  for  the  transverse 
loading  of<j(x,  y). 


Figure  5 

Truss-Core  Sandwich  Panel 
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The  reometry  of  the  structure  is  defined  by  the  following 
dimensions : 

b)  "a"  ar  "b"  are  the  overall  dimensions, 

b)  hr  is  the  thickness  of  the  core, 

c)  tj.  and  tc  are  the  thickness  of  the  face  plates  and  the  truas- 
core,  respectively, 

d)  0  is  the  angle  made  by  the  core  elements  with  the  vertical. 

In  a  large  number  of  problems,  the  specifications  will  call  for 
a  design  in  which  the  overall  dimensions,  "a"  and  "b",  have  been  pre¬ 
determined,  Hence,  it  is  not  proper  to  treat,  "a"  and  "b"  as  design 
parameters. 

Although  the  structure  is  shown  with  all  other  dimensions 
independent  of  location  in  the  plate,  this  nec-i  not  be  the  case.  If 
the  manner  in  which  the  transverse  loading  varies  with  position  can  be 
limited,  as  well  as  the  typo  of  supports  that  will  be  admitted  for  the 
plate  boundaries,  then  the  manner  in  which  the  interior  stress  field 
varies  will  also  be  limited.  In  such  cases  it  is  possible  to  introduce 
smaller  or  lighter  members  ir.  those  legioi.s  m  which  the  level  of 
stress  is  low  compared  to  the  same  type  members  which  are  in  regions 
of  higher  stress  levels.  For  this  work,  ail  core  members  are  assumed 
to  be  Identical;  all  face  plate  members  are  assumed  to  be  identical; 
and  hc  and  0  are  assumed  to  ;  o  constant  with  respect  to  any  position 
in  the  plate. 

In  this  investigation,  it  is  assumed  that  the  material  used  is 
homogeneous,  isotropic  and  linearly  elastic.  However,  the  facing 
material  and  the  truss-core  material  is  taken  to  be  different  materials. 
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Tha  mechanical  behavior  of  an  isotropic  linearly  elastic  solid  is 
governed  by  two  constants  and  ner.ee  two  material  design  parameters 
are  introduced,  in  addition  to  the  four  geometric  design  parameters 
hc,  tr,  tc,  and  0.  Note  that  for  panels  of  different  materials  for 
facing  and  core,  it  i.3  only  necessary  to  specify  the  known  ratio  of 
tho  material  properties  and  not  two  additional  parameters. 

It  might  be  well  to  emphasize  the  fact  that  any  optimum  design 
procedure  must,  undoubtedly,  include  orthotropic  material,  limited 
plastic  action  and  variation  of  some  of  the  geometric  parameters  with 
position.  In  this  light,  the  present  work  is  to  be  viewed  as  a 
first  step. 

C,  Failure  Criteria 

The  first  requirement  of  a  design  is  that  it  does  not  fail  on  the 
basis  of  the  following  criteria: 

a)  excessive  deflection 

b)  strength  deficiencies 

c)  instability 

Except  in  those  cases  for  which  there  are  very  strict  allowable 
deflection  tolerances,  an  excessive  deflection  failure  will  usually  be 
associated  with  some  amount  of  "plastic"  action  occurring  somewhere  in 
tha  structure.  To  design  for  an  excessive  deflection  failure,  there¬ 
fore,  requires  analysis  of  the  structure  in  those  cases  in  which  part  of 
the  structure  is  undergoing  plastic  action.  Since  the  latter  problem  is 
a  formidable  one,  a  much  more  restrictive  critorion,  completely  dis¬ 
allowing  any  plastic  action  or  allowing  only  some  arbitrarily  set  limit 
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of  plastic  action,  is  frequently  substituted  for  the  excessive 
defection  criterion. 

A  material  failure  may  be  said  to  occur  when  the  internal 
forces  exceed  the  strength  or  a  specified  stress  level  of  the 
material.  Once  again  the  material  will  usually  not  fail  without 
first  undergoing  some  plastic  action.  In  a  statically  indeterminate 
structure  like  the  truss  core  sandwich  plate  the  onset  of  plastic 
action  will  result  in  a  redistribution  of  internal  forces  which  will 
change  the  stress  distribution.  Once  again,  the  calculation  of  this 
redistribution  is  extremely  complicated  with  the  result  being  valid 
only  for  the  specific  loading  history  used  in  the  analysis. 

Two  other  material  failures  which  can  arise  but  are  difficult 
to  quantify  or  express  analytically  are  due  to  creep  and  fatigue. 
However,  these  failures  are  beyond  the  scope  of  this  work. 

An  instability  failure  occurs  for  a  given  loading  if  the 
distribution  of  internal  forces  calculated  for  that  loading  will 
undergo  a  significant  change  when  a  small  external  perturbing  factor 
is  introduced.  For  the  truss  core  sandwich  panel,  it  is  possible 
to  distinguish  between  a  local  instability  which  arises  due  to  an 
instability  in  the  response  of  an  individual  component  and  an  overall 
instability  which  arises  because  the  manner  in  which  the  components 
are  Joined  causes  the  articulated  structure  to  become  unstable.  It 
might  be  pointed  out  that  for  a  statically  indeterminant  structure 
a  local  instability  does  not  necessarily  mean  a  catastrophic  structural 
failure  but  may  merely  result  in  a  redistribution  of  internal  forces. 

For  the  truss  core  sandwich  panel  subjected  to  a  transverse 
loading,  there  is  no  possibility  of  an  overall  Instability  failure. 
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There  la,  of  course,  a  possibility  of  a  local  instability  occurring 
and  for  the  geometry  of  the  truss  core  treated  here,  such  instabilities 
will  be  catastrophic. 

One  typo  of  failure  that  sometimes  occurs  but  is  not  considered  in 
this  treatment  is  a  failure  of  the  Joints  of  the  structural  componente. 

In  fact  there  is  much  experimental  evidence  which  indicates  that 
improperly  designed  Joints  initiate  the  failure  of  a  structure,  however, 
a  detailed  consideration  of  this  problem  is  beyond  the  scope  of  this  work. 

For  the  truss  core  sandwich  panel  the  structure  will  be  said  to 
have  failedj  (a)  if  the  state  of  stress  anyvhero  in  the  structure  is 
beyond  the  linear  range  or  (b)  if  elastic  instability  occurs  in  any 
individual  component.  Thus  the  failure  modes  are  consistent  with  the 
analyses  developed  in  Chapter  1. 

D,  Failure  Modes 

Since  the  analysis  presented  in  Chapter  1  has  been  developed  for 
linear  elasticity  only,  the  expressions  for  describing  failure  modes 
will  have  to  be  consistent  with  that  development.  Further,  the  basic 
premise  in  the  development  of  the  load-displacement  relationships  for 
the  plate  elements  of  the  panel  has  been  the  one-dimensional  character¬ 
ization  (see  Appendix  1).  Consistent  with  this  approach,  then,  it  ia 
reasonable  to  define  failure  modes  in  terms  of  the  average  stresses 
IT',  p  ,  CT  ,  for  the  plate  elements. 

Due  to  the  nature  of  the  loading,  that  is,  laterally  applied  to  the 
truss  core  panel,  and  simply  supported  edges  it  is  obvious  that  the 
important  elements  to  consldor  an  I  t.hoir  primary  stresses  are  (see 
Figure  5):  (1)  upper  fa<  e  elements  sin  h  as  (lt-6)  subjected  to  biax.al 
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compression;  (2)  lover  face  elements  such  as  (5-7)  subjected  to  biaxial 
tension;  ana  (3)  web  elements  such  as  subjected  to  combined 

•trasses  of  compression,  shear  and  in-plane  flexure.  Since  in  any 
practical  case  it  is  not  always  obvious  what  the  degree  of  fixity  is  at 
the  edges  of  the  panel,  it  is  conservative  to  U6e  simply  supported  edges. 
Now  the  specific  criteria  for  these  elements  will  be  listed  and  as  such 
represent  the  constraints  which  are  Imposed  on  the  structure  and  which 
must  not  be  violated  when  determining  «  minimum  weight  design. 

The  upper  face  elements  are  subjected  to  essentially  biaxial 
compression  as  shown  in  Figure  6,  which  can  be  evaluated  as 


mot e  . 


Figure  6 

Stresses  on  Typical  Element 
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However,  it  can  be  shown  that  since  b  ,  ths  stability  of  this 

element  can  be  reduced  to  a  one  dimensional  problem.  Hence,  the 
buckling  criterion  is 


CT, 


AX 


Pil Ejr 
3(i-vf4) 


(2.3) 


so  that  the  constraint  introduced  for  the  stresses  in  the  upper  face  is 


t r  $ 


1_ 

3 


(s.b> 


The  lower  face  elements  are  subjected  essentially  to  biaxial 
tension  (Figure  6)  which  can  be  evaluated  as 


CTyy  =  ‘T^- 
ov.»  ” 

^  /? 


(2.5) 

(2.6) 


In  this  combined  stress  situation  it  is  usual  to  resort  some  theory  of 
failure  depending  on  the  type  of  materiel  used.  For  example,  if  the 
material  is  ductile  a  widely  accepted  thoory  is  Maximum  Shear  which 
interprets  failure  as  occuring  at  first  yielding  of  the  ductile  material; 
if  the  material  is  brittle  then  many  times  the  Maximum  Stress  Theory  is 
used.  Thus  the  constraint  introduced  for  the  stresses  in  the  lower  face 
is  simply  that  the  larger  principal  tensile  stress  is  less  than  or  equal 
to  the  fracture  stress  in  a  uniaxial  tensile  test  of  the  material. 

The  web  element  is  essentially  subjected  to  a  combination  of 
stresses  resulting  from  O'*.*.  ;  CTyy  and  CT*y  as  shown  in  Figure  6. 

In  terms  of  the  average  stresses 
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Cyy  =  — 

1 

2 

L 

te 

(2.8) 

Cr*y  =  - 

1 

2 

-fr 

(2.9) 

Equation  (2.7)  i»  the  O,,,  stress  developed  as  a  ccnseqv.ance  of  the  O' 
etresaes  and  ia  giTen  in  Appendix  1  by  Equation  (A. 7a).  Note  that 
since  the  analysis  given  in  Chapter  1  evaluates  £3  (*)  by  a  truncated 
aeriee  then  can  be  found  by  termwise  differentiation.  Although 

no  rigorous  theory  exists;  to  account  for  buckling  of  a  long  thin  plate 
under  coabined  loading.  Reference  $  recommends  a  criterion  of  the 
following  type : 


♦  R*  *  1  (2.10) 

applied  stress 

m  T"" 

allowable  stress 

1,2,3  index  which  indicates  the  type  of  loading,  e.g., 
compression'0),  shear'8),  or  flexure'") 

x,y,t  numerical  exponents. 

In  particular,  Reference  6  in  Figure  27  presents  interaction  curves 
for  various  combinations  of  compression,  bending  and  shear.  Thus,  a 
criterion  given  by  Equation  (2.10)  or  data  from  interaction  curves 
the  third  oonotraint  is  Introduced  which  cannot  be  violated  during  the 
optimisation  process. 
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5,  Optimisation  Procedure 

Based  on  the  fora  of  analysis  developed  in  Chapter  1,  it  is  readily 


apparent  that  for  ths  trust:  core  sandwich  panel  under  lateral  loads 
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optimization  by  extremizution  o'.  a  functional  form  Is  impractical  but 
rather  a  nume::? cnl  procedure  will  have  to  be  utilized.  Consequently, 
optimum  design  will  not  '^e  defined  as  that  one  in  which  the  most 
obvious  failure  modes  occur  simultaneously.  Instead,  the  criterion 
used  will  be  to  choose  those  valt.os  of  the  design  parameters  which 
.sake  the  total  vteight  a  minimum  subject  to  the  restriction  that  none 
the  failure  ifiteria  are  vi  ilated. 

7/ he  reccsiflie  tded  method  o*  achieving  an  optimum  design  will  be 
the  direct  approach  as  described  in  Section  A.  Since  no  analytical 
expressions  exist  at  this  time  which  predict  the  stress  distributions 
fof  the  panels  considered  in  this  report,  the  direct  approach  appears 
to  be  the  most  straight  forward  and  feasible  technique  and  in  essence 
concurs  with  the  philosophy  promulgated  by  R.  Bellman  in  Reference  7. 

The  two  material  parameters,  — £  and  <mci  the  four  geometrical 

Ec 

parameters,  hc ,  tc ,  tf,  and  0  ,  define  a  six  dimensional  space  with 
a  one  to  one  correspondence  between  a  point  in  space  and  a  particular 
design.  Since  every  point  In  the  space  obviously  will  not  correspond 
to  a  design  which  will  not  fail,  the  first  step  will  be  to  ascertain 
the  region  in  space  which  corresponds  to  debigns  which  do  not  fail. 
These  valid  designs  are  numerically  determined  by  the  procedure 
given  in  the  flow  chart,  Figure  7.  Thus,  the  hypersurface  separating 
the  region  of  no  failure  from  the  region  of  failure  is  obtained.  As 
a  consequence  of  these  calculations,  the  extremum  regions  on  the 
weight  hyperBurface  will  be  grossly  defined.  The  next  step  is  to 
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refine  the  geometrical  parameter  increments  tn  the  region  of  minimum 
panel  weight  u-.icil  a  design  *:oniJ  guration  is  secured  to  whatever 

accuracy  is  deemed  necessary. 

Referring  to  Figure  5,  the  entii  weight  of  the  structure 

(q^cluBive  of  weight  of  adhesives  or  fas  eners),  expressed  in  terms  of 

tye  deafen  parameters,  is  easily  seen  to  t3 

+  2-Pftf]  <2 

W^ero  and  are  the  weight  densities  of  the  core  and  the  faces, 
respectively,  and  all  other  quantities  are  defined  in  Figure  5.  It  is 
immediately  seen  that  neither  nor  the  ratio  of  mechanical  properties 
enters  directly  into  the  weight,  however,  it  would  be  erroneous  to 
conclude  from  this  that  the  depth  d--'-  not  enter  into  the  weight. 

Only  two  of  the  three  parameters,  b,  h  ,  and  9  can  be  chosen  to  be 
independent  for  the  geometry  shown  and  Equation  (2.12)  takes  b  and  9 
as  the  independent  parameters.  Th8  ratio  of  the  mechanical  properties 
will  also  enter  the  picture  in  an  indirect  manner  since  they  will  have 
an  effect  on  what  values  for  the  other  design  parameters  constitute 
a  safe  design. 
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APPE'DIX  1 

GOVEPNI?£»  SQUAT;  -‘g  ft*  QN?,-PIHRNS10NAL  ELEMENTS 

It  is  the  purpose  of  this  appendix  to  investigate  the  validity 
of  p  °*ulating  that  the  basic  elements  have  only  one  significant 
dimension  and  to  obtain  '!,iations  governing  the  response  of  these 

one  dimensional  models. 

The  restriction  of  the  basic  elements  to  thin  plates  (i.e,  vith 
respect  to  the  other  two  dimensions  and  the  rate  of  spatial  variation 
of  the  loading)  subjected  to  in-plane  loading  allows  average  values  to 
be  used  a a  far  aa  the  thickness  direction  is  concerned.  The  resulting 
two  dimensional  theory,  generalised  plane  stress  theory,  is  well  under¬ 
stood  and  the  requirements  for  its  validity  need  not  bo  investigated 
here.  It  is  the  further  reduction  from  a  two-dimensional  problem  to  a 
one-dimensional  problem  that  needs  to  be  investigated. 

Case  I: 

Consider  tho  boundary  value  problem  shown  in  Figure  Al. 


Figure  Al 

with  the  boundary  conditions  at  x  •  0  and  x  -  1  not  being  specif 


uJ 


as  yet. 
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equations  in  conjunction  with  the  bou .  conditions  as  listed  in 
(A. U)  is  not,  in  general,  possible.  however  if  "h1-  is  much  less  than 
any  characteristic  dimension  in  the  x  direction  (i.e.  either  the  length 
or  some  measure  of  the  rate  of  variation  of  <rCn){ir  pCx>);  then  an 
approximate  solution  can  be  achieved  on  the  basis  of  the  following 
assumption:  the  y  variation  of  the  unknown  displacements  U.Cx,y) 
and  )i»  expandable  in  a  power  series.  Under  the  assumption  of 

small  h  relative  tc  dimensions  in  the  x  direction  these  power  series 
expansions  will  converge  quite  rapidly  hence  a  good  approximation  can 
be  obtained  by  truncating  after  the  first  few  terms.  Realising  this, 
the  truncation  is  introduced  in  the  very  beginning  rather  then 
obtaining  first  the  exact  solution  in  series  form  and  then  truncating. 
It  is  important  to  notice,  however,  that  the  truncated  series  is 
only  approximate  and  as  such  cannot  satisfy  the  pair  of  partial 
differential  equations  discussed  above  but  rather  can  satisfy  only 
certain  aspects  of  them. 

The  form  of  the  loading  applied  to  the  edges  at  y  *  t  h/2  will 
give  rise  to  a  u(*,y)  which  is  antisymmetric  in  y  and  a  V  (x,y) 
which  is  symmetric  in  y.  Noting  this,  thon  if  the  power  series 
expressions  is  truncated  after  one  tern  the  following  forn  for  the 
displacement  field  results 
a  U.y')  = 
ir  (x  ,y  )  *  Y (X) 

where  Q (x)  and  If  GO  are  to  be  determined  by  satisfying  some  aspects 
of  the  governing  equation.  The  strain  ann  stress  fields  associated 
with  the  above  displacement  field  is  obtained  by  direct  substitution 
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In  addition  to  satisfying  the  equilibrium  equations  in  an 
avers gs  sense,  the  requirement  that  stress  couple  equals  sero  can 
also  be  satisfied.  Multiplying  each  of  Equation  (A.l)  by  y  and  then 
integrating,  again  making  use  of  Equation  (A.ii)  the  results  arei 


(A. 10  a) 


(A. 10  b) 


By  direct  substitution  of  Equation  (A. 7)  into  Equation  (A. 10),  it  is 
readily  seen  that  Equation  (A. 10  b)  is  identically  satisfied  while 
Equation  (A. 10  a)  introduces  the  following  restriction  on  {3  {*) 
and  T  (*)  t 


ZGt 


L  h 


-  Jh_ 

3(l-v)dxt 


(A. 11) 


Equations  {A, 9)  and  (A.  11)  f»tvr  the  deMred  relations  between 
the  loading  and  the  displacement  field  subject  to  the  restriction 
already  noted,  A  sore  rigorous  analysis  using  energy  concepts  would 
show  that  It  is  logically  consistent  to  introduce  a  shape  factor 
(terawd  hj'1)  into  Equation  (A. 7  c).  That  is, 


cr*.  y 


dX\ 
d*  ) 


(A, 7  d) 


Introducing  this  factor  permits  Equations  (A. 7  d)  and  (A, 11)  to  be 
written  as 


5) 


n\)<  -Aair-nu 


(X (x)  =  -2  K*6  t  -  K2G  h £  dL&ti. 

d  x  c/a1 

pU)«[4KlGt  _2G hi.  dl.l|3(;x')4.2.KaGtdlL>0A*l2) 
'  *-  h  3(/-i?)dx*J'  dx 


A  word  should  be  said  about  the  boundary  conditions  to  be  applied 
to  the  ends  s  »  0  and  x  ■  J.  Obviously,  it  is  not  possible  to  satisfy 
the  exact  boundary  conditions  of  the  three  dimensional  theory  but  only 
satisfy  these  boundary  conditions  in  an  average  sense.  9 y  virtue  of 
St,  Venant's  principle,  however,  it  la  known  that  the  error  associated 
with  not  satisfying  theeo  conditions  exactly  is  confined  to  the  region 
in  the  vicinity  of  the  two  ends. 


Case  lit  Next,  consider  the  boundary  value  problem  shown  in 
Figure  A2, 

_  (y _ \t  ?<*) _ 


Figure  A 2 


with  the  boundary  conditions  at  x  •  0  and  x  •  Ji  being  left 
unspecified  for  the  tins  being. 

The  governing  equations  are  again  (A.l)  through  (A. 3)  with  the 
boundary  condition*  on  y  •  t  h/2  being  expressed  aatheoaticeUy  as 

crTy  = 
ry  zt 


(a. u) 


crv*  a  o 
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A  solution  of  the  following  form  is  taken 

L(  (x,y)  --  U  '  (x,y ) 

I/' (x.y)  -  v;  My  +  ir'(x,y) 


where  it  has  been  assumed  that 


U7(x)  -  0r^)rcx) 
4Gf 


(a. 15) 


By  substituting  this  form  of  the  solution  in  Equations  (A.l)  through 
(A. 3)  *nd  (A. 13)  results  in  the  following  boundary  value  problem  on 
'<X  y )  and  cr  *6r,  y  )  . 


Equilibrium : 


cr  +•  cr*.sy,  y 

O^y,  A  +  OVy>v 


f  -  ^  dx 

Jj jj2)  _d\£  ooy 
(\-u)  c/x1 


(A. 16) 
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Boundary  Conditions  on  y  *  I  h/?i 
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(A. 17) 


(A. 18) 


(A. 19) 
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where  the  sinus  sign  is  cor  roc  t  for  y  -  4  h/2  and  the*  plus  sign 
i*  correct  for  y  «  -  h/2. 

Mow,  without  actually  solving  tha  boundary  value  problem  for 
U.'(Kty)  am}  yj  ,  it  is  desirable  to  obtain  soma  knowledge  of  the 

order  of  magnitude  of  these  terms  relative  to  dZ  (x)  y 
To  do  this  properly  all  dimensions  are  first  normalised  with  respect 
to  some  characteristic  length,  say  h. 


Let 


(A. 20) 


In  terms  of  and  ^  as  the  independent  variables  the  boundary  value 
problem  on  U.’Cf »nd  4r'(£(  hy)  becomes 
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In  Equations  (A.?0)  through  (A ,2h) ,  th*>  double  prime  has  been  intro¬ 
duced  to  emphasize  the  fact  that  the  independent  variables  have  been 
changed* 

It  can  be  shown  that  all  of  the  forcing  terms  on 
and  [V"(£f  »\)A]  are  dependent  on  darivate9  of  if, ' 

with  respect  to  £  .  Since  this  is  true  and  since  the  first  portion 

of  the  solution  is  equal  to  U£(£)J^i^it  is  consistent  to  neglect  the 
second  portion  relative  to  the  first  provided 


<<  1 


vUS) 

The  result  of  doing  this  is  to  obtain  sb  a  solution  the 


following 
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which  is  used  in  Chapter  1. 

Nothing  was  said  of  the  boundary  conditions  applied  at  *  *»  o 
and  x  *  Jl  ,  Once  again  under  the  assumption  of  h  <*■&  ,  the  effect 
of  the  boundary  will  be  small  except  near  the  extremities  of  the  member. 


Case  lilt  The  final  case  is  the  boundary  value  problem  ae  shown  in 
Figure  A3. 


~  rrfx) 

Figure  A3 


with  the  conditions  at  v.  •  0  and  x  •  1  left  unspecified  for  the  noaent. 
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A  a  before  it  can  be  shown  that  If  h  <<  M.  and  much  less 
than  any  rate  of  spatial  variation  of  '??''<*)  then  it  is  justifiable  to 
deal  with  an  average  value  relative  to  the  y  direction.  Such  an 
attack  will  result  in  the  number  being  treated  as  one  dimensional. 
The  solution  is 
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